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Abstract Multiple recent investigations of solar magnetic field measurements 
i have raised claims that the scale-free (fractal) or multiscale (multifractal) pa- 

rameters inferred from the studied magnctograms may help assess the eruptive 
potential of solar active regions, or may even help predict major flaring activity 
stemming from these regions. We investigate these claims here, by testing three 
f>r^ I widely used scale-free and multiscale parameters, namely, the fractal dimen- 

sion, the multifractal structure function and its inertial-range exponent, and the 
. turbulent power spectrum and its power-law index, on a comprehensive data 

I set of 370 timeseries of active-region magnctograms (17, 733 magnctograms in 

total) observed by SOHO's Michelson Doppler Imager (MDI) over the entire 
Q^, Solar Cycle 23. We find that both fiaring and non-flaring active regions exhibit 

^ ' significant fractality, multifractality, and non-Kolmogorov turbulence but none 

' of the three tested parameters manages to distinguish active regions with major 

C/5 , flares from flarc-quict ones. We also find that the multiscale parameters, but 

^ ' not the scale-free fractal dimension, depend sensitively on the spatial resolu- 

tion and perhaps the observational characteristics of the studied magnctograms. 
^ , Extending previous works, we attribute the flare-forecasting inability of fractal 

' and multifractal parameters to i) a, widespread multiscale complexity caused by 

. a possible underlying self-organization in turbulent solar magnetic structures, 

flaring and non-flaring alike, and ii) a. lack of correlation between the fractal 
properties of the photosphere and overlying layers, where solar eruptions oc- 
cur. However useful for understanding solar magnetism, therefore, scale-free and 
multiscale measures may not be optimal tools for active-region characterization 
in terms of eruptive ability or, ultimately, for major solar-flare prediction. 
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1. Introduction 



The ever-increasing remote-sensing capabilities of modern solar magnetographs 
have led to the undisputed conclusion that solar (active region in particular) 
magnetic fields exhibit an intrinsic complexity. "Complexity" is a term commonly 
used to describe an array of properties with one underlying characteristic: a 
scale-invariant, self-similar (fractal) or multiscale (multifractal) behavior. The 
measured photospheric magnetic fields in active regions are indeed multifrac- 
tal {e.g. Lawrence, Ruzmaikin, and Cadavid 119931 lAbramenkoi 12005^ . that is, 
consisting of a number of fractal subsets. As such, they are also fractal, with a 
fractal dimension equal to the maximum fractal dimension of the ensemble of 
fractal subsets. 

Fractality is a mathematical property but with important physical impli- 
cations. Scale-free or multiscale manifestations are thought to stem from an 
underlying self-organized, or self-organized critical (SOC), evolution in active 
regions. Self-organization refers to the internal, intrinsic reduction of the vari- 
ous parameters (also called degrees of freedom) of a nonlinear dynamical sys- 
tem, such as a solar active region, into a small number of important parame- 
ters that govern the system's evolution and, perhaps, its dynamical response 
( [Nicolis and Prigogine, 1989[ ). Assumptions on the nature of just these impor- 
tant parameters can lead to models of active-region emergence and evolution 
encapsulated in simplified cellular automata models ( Wentzel and Seiden, 1992} 



Seiden and Wentzel, 1996||Vlahos et al, 20021 [Fragos, Rantsiou, and Vlahos, 2004D .| 



Self-organized criticality, on the other hand, implies that the self-organized sys- 
tem evolves through a sequence of metastablc states into a state of marginal sta- 
bility with respect to a critical threshold. Local excess of the threshold gives rise 
to spontaneous, intermittent instabilities lacking a characteristic size ( |Bak, Tang, and Wiesenfeld] 
[II87; Bak, 1996,) . 

The intrinsic self-organization in solar active regions may be attributed to 
the turbulence dominating the emergence and evolution of solar magnetic fields. 
Tangled, fibril magnetic fields rising from the convection zone can be explained 
via Kolmogorov's theory of fluid turbulence {e.g., [Brandenburg et al\ 119901 



[Longcope, Fisher, and Pevtsov[ 119981 [Cattaneo, Emonet, and Weiss 120031 and 
others). Turbulence in the generation and ascension of solar magnetic fields 



leads to turbulent photospheric fiows {e.g., Hurlburt, Brummel, and Toomre 
I1995p . Thus, the turbulent photosphere is viewed as a driver that gradually but 
constantly perturbs an emerged magnetic-fiux system, such as an active region, 
dictating self-organization in it and possibly forcing it toward a marginally sta- 
ble, SOC state {e.g., Vlahos and Georgoulis| [2004)) . Turbulent action does not 
cease in the photosphere, but it extends into the solar corona. However, coronal 
low-/3 turbulence may not be the Kolmogorov fluid turbulence applying to the 
high-/3 plasma of the convection zone and the photosphere. Instead, it might 
be an intermittent magnetohydrodynamic (MHD) turbulence ( [Kraichnan, 1965} 



Biskamp and Welter, 1989). 



Fractal, multifractal, and turbulent properties of photospheric active-region 
magnetic fields have been intensely studied in recent years. Fractality is tradi- 
tionally investigated via the fractal dimension, often inferred using box-counting 



SOLA_ms.tex; 4 January 2011; 1:54; p. 2 



Are Flaring Active Regions More Complex than Others? 



3 



techniques {e.g., IMandelbroti I1983p . Box-counting is also used for multifrac- 
tal studies in space and time {e.g., lEvertsz and Mandelbrotl I1992p . involving 
also generalized correlation dimensions ( [Georgoulis, Kluiving, and Vlahos, 1995] 



Kluiving and Pasmanter, 1996). A commonly used multifractal method that does| 



not require box counting is the calculation of the multifractal structure function 
spectrum ( [Frisch, 1995[ ). Moreover, a practical method for quantifying turbulence 
is the calculation of the turbulent power spectrum, stemming from the original 
work of |Kolmogorov| (|194ip . If the power spectrum shows a power law over a 
range of scales, perceived as the turbulent inertial range, its slope determines 
whether the inferred turbulence is Kolmogorov-like (scaling index w —5/3) or 
Kraichnan-like (scaling index w —3/2) if either of these two applies. 

Multiple studies on fractality, multifractality, and turbulence in photosphcric 
active-region magnetic fields have raised claims that flaring active regions exhibit 
distinct, distinguishable complexity. These works might lead to the impression 
that fractal, multifractal, or turbulent measures hold significant flare-predictive 
capability or, at least, they might be used to identify flaring active regions 
before they actually flare. To summarize some of these works. Rbramenko et al.\ 
(|2003p suggested that a "peak in the correlation length might be a trace of 
an avalanche of coronal reconnection events" . [McAteer, Gallagher, and Ireland] 
(|2005p reported that "solar flare productivity exhibits an increase in both the 
frequency and GOES X-ray magnitude of flares from [active] regions with higher 
fractal dimension" . Further, lAbramenkol (I2005P found that "the magnitude of 
the power index at the stage of emergence of an active region ... reflects its future 
flare productivity when the magnetic configuration becomes well evolved" , while 
GeorgouUsI (j2005p reported that "the temporal evolution of the [inertial-range] 



scaling exponents in flaring active regions probably shows a distinct behavior 
a few hours prior to a flare" . More recently, IConlon et al. \ (I2008P worked on a 
sample of four active regions and reported evidence for a "direct relationship 
between the multifractal properties of the flaring regions and their flaring rate" , 
while IHewett et al.\ (|2008p . reporting on "preliminary evidence of an inverse 
cascade in active region NOAA 10488" found a "potential relationship between 
energy [power-spectrum] scaling and flare productivity" . Many of these works 
are also reviewed by |McAteer, Gallagher, and Conloii] (|2010p . 

If the above findings are confirmed, they may well lead to notable improve- 
ments in our physical understanding of active regions and in highlighting possible 
differences between flaring (that is, hosting major flares) and non-flaring (that 
is, hosting only sub-flares) regions. In IGeorgouiisl ([2005P we studied three dif- 
ferent scale-free and multiscale parameters, namely, the fractal dimension, the 
spectrum of generalized correlation dimensions, and the structure-function spec- 
trum and its inertial-range exponents over a limited magnetogram sample of six 
active regions, three of them hosting at least one major flare (M- or X-class 
in the GOES X-ray 1-8 A flare classiflcation scheme). In one case of a X- 
flaring active region - NOAA active region (AR) 10030 with an X3 flare at the 
time of the observations - we noticed a sharp preflare increase of the inertial- 
range exponent of the structure functions followed by a significant (« 20% and 
much above uncertainties), permanent decrease after the flare. Wc suggested 
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that this analysis should be repeated on a much larger sample of both flaring 
and non-flaring regions to determine whether this behavior was incidental. 

In this study we analyze a comprehensive sample of 370 timeseries of active- 
region magnetograms, with each timeseries corresponding to a different active 
region. In this sample, 77 active regions hosted at least one M- or X-class flare 
during the observations and they are considered flaring (17 X-class flaring, 60 
M-class flaring), while the remaining 293 active regions were not linked to major 
flares and are hence considered non-flaring. We calculate three of the most 
promising scale-free and multiscale measures on this data set, namely, the frac- 
tal dimension, the multifractal structure function spectrum, and the turbulent 
power spectrum. A detailed description of the data and techniques follows in 
Section [2] In Section [3] we test the sensitivity of the calculated parameter values 
on the spatial resolution of the studied magnetogram. A statistical analysis of 
the active-region sample is performed in Section 2] while Section [5] summarizes 
the study, discusses the results, and outlines our conclusions. 



2. Design of the Study 
2.1. Magnetogram Data 

Our active-region sample has been constructed using data from the Michelson- 
Doppler Imager (MDI: [Scherrer et aZ.[[TD95p . onboard the Solar and Heliospheric 
Observatory (SOHO) mission. We acquired the entire MDI magnetogram archive 
from mid-1996 to late-2005. The archive consists of full-disk line-of-sight solar 
magnetograms taken at a 96-minute cadence with a linear pixel size of « 1.98 
arcsec (a mean « 1440 km at solar disk center, depending on Sun-Earth dis- 
tance). This analysis uses purely line-of-sight. Level 1.5 SOHO/MDI magnetic 
field measurements that are known ( Berger and Lites, 2003[ ) to underestimate 



sunspot and plage fields (more recent. Level 1.8.2 sensitivity corrections to the 
MDI full-disk magnetograms are not used in this study because our full-disk mag- 
netogram dataset was constructed in 2007 and the recalibrated magnetograms 
were posted in December 2008). Nonetheless, we have avoided applying any 
additional corrections to the data to avoid a possible impact on the morpholog- 
ical characteristics of the regions studied, since fractal and multifractal analysis 
highlights exactly these characteristics. To reduce the impact of projection effects 
acting on the magnetic field vector, we restrict our study to a 60° longitudinal 
region centered on the central solar meridian. Use of this zone introduces a 
systematic underestimation in the normal magnetic field component by a factor 
up to w (l-cos(0)) ~ 0.14, or 14%, at the Equator, for a central meridian distance 
9 = 30° approximated by the angular difference between the local normal and 
the line of sight for an observer at Earth. Conventionally, this underestimation 
factor is deemed tolerable when the line-of-sight field component is used as 
a proxy of the normal field component. Within this 60° meridional zone we 
identified and extracted active regions using our automatic active-region iden- 
tification algorithm (ARIA) , detailed in |LaBonte, Georgoulis, and Rust| (I2007P 
and in [Georgoulis, Raouafi, and Henney| (j2008p . Our ARIA extracts portions of 
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Figure 1. Pictorial output of our ARIA, applied to a full-disk SOHO/MDI magnctogram 
acquired on 29 October 2003. A 60° meridional zone centered on the central meridian is 
indicated by the thick dashed brackets. Three active regions fulfilled our selection criteria, 
namely, NOAA ARs 10486, 10487, and 10488. The portion of the disk found to include each 
region is shown by the thin dashed circles - the actual portion extracted for each region is shown 
by the circumscribed squares. The NOAA labels for each region are provided automatically. 
The flux-weighted centroids for each active region are represented by white crosses. 

the solar disk corresponding to active regions by means of pattern recognition 
in which the unit length is one supergranular diameter (40 arcsec). A typical 
example is shown in Figure [TJ An active region is chosen for further study if 
its flux-weighted centroid, shown by the white crosses in the selected regions 
of Figure [1] falls within the above-mentioned 60° meridional zone. Notice, for 
example, that NOAA AR 10487 is selected in Figured] because its flux- weighted 
centroid lies within the above zone; parts of it, however, extend beyond this area. 

Besides the automatic active-region selection process, each selected magnc- 
togram (out of a total of 17,733) was manually examined to exclude portions 
of other active regions that might intrude in the field of view. For example, 
NOAA AR 10486 in Figure [1] is included in its selection circle, but the square 
circumscribed on this circle crops sizable parts of NOAA ARs 10489 and 10491 
in its northwestern edge. These parts have been excluded in the subsequent 
analysis. Generally our ARIA performs quite well in distinguishing active re- 
gions but few incidences such as the above have been noted, especially in cases 
of densely populated active-region complexes, or "nests" , such as the one be- 
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Figure 2. Temporal distribution of our sample of 370 active regions over Solar Cycle 23. 
Shown arc the median observation times (color symbols) of each region together with the 
monthly-averaged sunspot number (curve). 



lieved to have occurred during the October-November "HaUoween" 2003 period 
( |Zhou et ai, 2007| . For our analysis, ARIA uses a maximum tolerated magnetic- 
flux imbalance of 50% in a given active region and a minimum active-region linear 
size of one supergranular diameter. For each of the 370 selected active regions 
we created a timeseries consisting of up to « 60 magnetograms taken every 96 
minutes corresponding to the approximately four-day period needed for each 
active region to traverse the 60°-analysis zone. 

To document the major flare history for each active region we browsed i) 
NOAA's GOES X-ray archive and ii) the Yohkoh/HXT flare catalog (available 



online, at http://gedas22.stelab.nagoya-u.ac.jp/HXT/catalogue/|. From the total 



of 370 active regions, 77 were unambiguously found to have hosted at least one 
M-class or X-class flare while within ±30° of the central meridian, with a total of 
24 X-class flares and 87 M-class flares. Our active-region sample roughly covers 
Solar Cycle 23, as shown in Figure [5] The solar cycle is represented by a 5-point 
running mean of the monthly-averaged sunspot number obtained by the Solar 
Influences Data Analysis Center (SIDC) of the Royal Observatory of Belgium. 

In addition to the above SOHO/MDI sample, our analysis includes three 
nearly simultaneous magnetograms of NOAA AR 10930, observed on 11 Decem- 
ber 2006. The line-of-sight components of these magnetograms are depicted in 
Figure [31 Figure shows the Level ID magnetogram (preferred over Level 2 
data in order to better qualify for comparison with SOHO/MDI Level 1.5 data) 
acquired by the Spectropolarimeter (SP: |Lites, Elmore, and Streander" I200ip of 
the Solar Optical Telescope (SOT) onboard the Hinode satellite and has a very 
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Figure 3. Nearly simultaneous, coaligned magnet ograms of NOAA AR 10930, acquired on 11 
December 2006: observations arc from (a) Hinode's SOT/SP, (b) SOHO/MDI high-resolution, 
partial disk magnetograph, and (c) SOHO/MDI full-disk magnctograph. The linear pixel sizes 
are 0.158 arcsec, 0.605 arcsec, and 1.98 arcsec for (a), (b), and (c), respectively. Shown is the 
line-of-sight magnetic field component saturated at ±2.5 kG (a), ±1 kG (b), and ±1.8 kG (c). 
Tic mark separation in all images is 10 arcsec. 



high spatial resolution 0.32 arcsec) with a linear pixel size of w 0.158 arcsec. 
Figure [3Jd shows the respective magnetogram taken by the SOHO/MDI high- 
resolution, partial-disk magnetograph, with a coarser linear pixel size of 0.605 
arcsec. Figure |3J: shows the Level 1.5 SOHO/MDI magnetogram extracted by 
a full-disk measurement with a much coarser linear pixel size of 1.98 arcsec. 
The three magnetograms have been initially coaligned by means of the pointing 
information provided separately for each. To further correct and deal with small 
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pointing inconsistencies, coalignment has been completed by a rigid displacement 
(translation) over the E-W and the N-S axes. Displacements are determined by 
the peak of the cross-correlation function between a given pair of images, with 
cross-correlation functions inferred by means of fast Fourier transforms. 

The three distinctly different spatial resolutions of the magnetograms of Fig- 
ure [3] will be useful when testing the sensitivity of our scale-free and multiscale 
parameters to varying spatial resolution (Section [3]). 

2.2. Scale-free and Multiscale Techniques 

The first parameter that we calculate is the scale-free, two-dimensional fractal 
dimension [Dq] of the active-region magnetograms. To calculate Dq we cover 
the magnetogram ficld-of-view with a rectangular grid consisting of square boxes 
with linear size [A] and area A x A. Assuming that the field of view is a square 
with linear size L and area L x L, the number of boxes needed to cover it is 
(i/A)^. Each of the boxes will have a dimensionless area e^, where e = \/L. 
Of the total (i/A)^ boxes we count those that include part of the boundary 
of a strong-field magnetic configuration (see below for the adopted strong-field 
definition). Then, varying the box size A or, equivalently, the dimensionless size 
e, we obtain different numbers [iV(e)] of information-carrying boxes. Correlating 
the various numbers N{e) with the respective box sizes e, we obtain the scaling 
relation 

iV(e) oc (l/£)^° . (1) 

For a non-fractal, Euclidean structure embedded on a plane we have N{e) — 
(L/X)^ ~ (1/e)^, so Dq = 2. If Dq < 2, we have fractal structures with a scale- 
free, incomplete filling of the field of view. The stronger the departure of Dq 
from its Euclidean value of 2, the finer the structure exhibited by the studied 
magnetic configuration. For Dq < 1 in a two-dimensional fractal, the structures 
are typically scattered into a scale-free hierarchy of small "islands', resembling 
what is known as fractal dust (e.g. ISchroederl |199ip . 

We infer the fractal dimension Dq by a least-squares best fit of the scaling 
relation of Equation ([T]). The uncertainty associated with the value of Dq is 
equal to the uncertainty of the regression fit. To guarantee a reliable inference 
of Dq, we demand that the dynamical range represented by the least-squares fit 
exceeds one order of magnitude. A very small fraction of magnetograms of non- 
fiaring active regions (~ 0.6%, or 111 magnetograms) happen not to comply with 
this requirement because of the regions' simplicity and scattered configurations; 
these magnetograms have been excluded from the analysis. Nonetheless, each of 
the 370 active regions of our sample fulfills the requirement with at least one 
magnetogram. 

|McAteer, Gallagher, and Ireland] (|2005p . relying on a substantial data set of 
w 10^ active regions, first reported that fiaring regions have fractal dimensions 
Dq > 1.2. Their finding was statistical, of course, meaning that the Dq > 1.2 
condition should be viewed as a necessary, but not sufficient, condition for ma- 
jor fiarc productivity. They also concluded, and we test their result here, that 
intensely fiaring active regions showed statistically higher fractal dimensions. 
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To better compare with the results of |McAteer, Gallagher, and Ireland] (|2005p . 
we follow their criterion when outlining the boundaries of active regions: first, 
we use a threshold of 50 G in the strength of the line-of-sight field component 
in order to define the outer contours of strong-field magnetic patches. Then we 
impose a lower limit of 20 pixels for the length of each contour, thus rejecting 
very small patches that could as easily belong to the quiet Sun. As an additional 
condition, we impose a lower flux limit of 10^° Mx for each patch. For the MDI 
low-resolution data this is nearly equivalent to saying that at least 100 pixels 
within the patch should have a line-of-sight field strength of at least 50 G, which 
is our threshold. Of course, a selected patch can contain fewer than 100 (but 
more than 20) pixels but with larger field strength in order to satisfy the flux 
condition. 

The second parameter that we calculate is the multifractal structure func- 
tion spectrum ( [Frisch, 1995[ ). The spectrum is given by 

5,(r) = (|$(x + r)-$(x)|«) (2) 

and it does not rely on box-counting or thresholding, contrary to Dq. Instead, on 
the magnetic- flux distribution [^(x)] of the active-region photosphere we define 
a displacement vector [r], also called the separation vector, and calculate the 
variation of the flux at this displacement. The variation is then raised to the 
power q, where g is a real, preferably positive number called the selector. Spatial 
averaging (()) of the structure function over x and all possible orientations of r 
gives rise to a unique, positive value Sq{r) of the structure function for a given 
pair (r = |r|, g). The resulting spectrum involves a range of r- values and a fixed 
value of q; different spectra are obtained for different g-values. 

The multifractal structure function is designed to highlight the intermit- 
tency present in a magnetic- flux distribution ([Abramenko et all 120021 l2003p . In 
the case of a multifractal, intermittent flux distribution, the structure function 
[Sq{r)] exhibits a power law 

Sqir) cx r«(9) (3) 

within a range of displacements, often referred to as the turbulent inertial range. 
The upper and lower extremes of the r-range correspond to, respectively, the 
maximum size of structures entering the inertial range and the scale over which 
ideal cascading of energy to smaller scales breaks down by dissipative effects. 
Higher values of the inertial-range scaling index ^(q) indicate a higher degree of 
intermittency, with (^{q) = q/i implying absence of intermittency. 

[Abramenko, Yurchyshyn, and Wang|(|2008p studied the structure function spec-| 
trum and implemented an additional suggestion by IFrischI (|1995p to examine 
the flatness function / = 56(r)/5'3(r)^ cx r"'', where 5 is viewed as the inter- 
mittency index with higher values implying a higher degree of intermittency. 
They concluded that photospheric and coronal magnetic fields are both in- 
termittent, with the intermittency in the photosphere preceding that in the 
corona and the corona responding to photospheric increases of intermittency. In 
their analysis, they used data from the same active region we show in Figure |31 
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namely, NOAA AR 10930, observed by ffinode/SOT/SP and SOHO/MDI high- 
resolution magnetographs, although not on the same day with the data used in 
this study. 

|GeorgoriIis| (|2005p used the structure function Sq (r) of Equation and the 
inertial-range scaling exponent to show that the photospheric magnetic 
fields of the six active regions of the study were indeed multifractal and intermit- 



tent, departing strongly from ({q) = q/3. More importantly, Georgoulis (12005^ 
presented an example of a X3 flare that occurred in NOAA AR 10030, where 
the intermittcncy peaked « 1 — 2 hours prior to the flare and decreased sharply 
after the flare. The decrease was permanent and C(<z) continued decreasing w 1.5 
hours later, when an Ml. 8 flare also occurred in the active region (Figure 7 of 
IGeorgorillsl I2005P . The change in the degree of intermittcncy was best seen for 
a selector q G (3,3.5). Here we investigate whether this distinct behavior, seen 
in only one example, is part of a systematic tendency. To this purpose we study 
the temporal evolution of the scaling index ({q — 3), inferred by a least-squares 
best fit of the scaling relation of Equation ([3]) . The uncertainty of the regression 
fit is treated as the uncertainty of the value of C(3). For the uncertainty of the 
C(3)-differcnccs we simply propagate the uncertainties of the two C(3)-values 
that make these differences. 

The third parameter we study is the turbulent power spectrum [E{k)]. In 
the case of a turbulent fiux system, there exists an incrtial range of wavcnumbers 
k reflected on a power-law form for E{k): 

E{k) oc fc-" , (4) 

where a is the inertial-range exponent. We demand at least one order of mag- 
nitude as the dynamical range of the least-squares best fit used to infer a and 
we attribute the uncertainty of the fit to the uncertainty of a. Like Sq{r) and 
unlike Dq, no box-counting or thresholding is required to infer a. 

lAbramenkol (|2005p suggested that the fc-range corresponding to length scales 
r e (3, 10) Mm should be scaled according to the representative a-value for 
high-resolution MDI magnetograms. Although we use MDI low-resolution data 
for this analysis, we follow lAbramenkol s suggestion and make sure that the 
(3, 10) Mm scale range is included in the power-law fit. We also attempt to test 
the conclusion of lAbramenkol (|2005p . relying on a sample of 16 high-resolution 
MDI active-region magnetograms, that the inertial-range exponent a rcficcts 
the future flare productivity of an active region, with larger a-valucs implying a 
higher flare index (Figure 8 of lAbramenkol I2005P . It is also important to mention 
here lAbramenkol s (j2Q05p conclusion that the scaling index a is not particularly 
useful for the prediction of imminent flares in the studied active regions. 

A typical calculation of the fractal dimension Dq, the inertial-range scaling 
index and the scaling exponent a of the turbulent power spectrum for the 
low-resolution MDI magnetogram of Figure |3t is depicted in Figure HI 

With the exception of the turbulent power spectrum, |GeorgouIis| (|2005p used 
all of the above methods, including the spectrum of multifractal generalized 
correlation dimensions. The latter showed that the six active regions of the 
sample had clearly multifractal photospheric magnetic flelds. The method is not 
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Figure 4. Typical calculation of the scale-free and multiscalc parameters that will be used 
to quantify the complexity of solar active regions. The example refers to the low-resolution 
SOHO/MDI magnetogram of Figure [3};, shown in (a). The magnetogram is saturated at ±1 
kG. Contours indicate the areas where the linc-of-sight field strength exceeds the threshold of 
50 G. Tic mark separation is 10 arcsec. Also shown are (b) the scaling relation of Equation ((TJ 
yielding the fractal dimension Do, (c) the structure function Sq—3{r) of Equation ((JJl, yielding 
the inertial-range scaling index C{3), and (d) the turbulent power spectrum of Equation (|4j, 
yielding the inertial-range scaling exponent a. The two dotted lines indicate the wavenumbers 
of the desired length scales of 10 Mm (fc < 1) and 3 Mm (fc > 1). Gray curves in b, c, and 
d correspond to N{e), 5g=3(r), and E{k), respectively, while the blue lines are the respective 
least-squares best fits. 



used in this study, however, because in |GeorgoiiIIs| (|2005p it clearly failed to 
distinguish between flaring and non-flaring regions even for the handful of active 
regions studied. 

The objective of this work is to test the flare-predictive capability of the 
above - reported to be promising - parameters. This capability will be tested by 
comparison with the capability of a standard, traditional parameter reflecting 
the size of active regions, namely, the unsigned magnetic flux 

*tot = / |S„|d5 , (5) 
Js 

where Bn is the normal magnetic field component (approximated by the line- 
of-sight field component near disk center) over the magnetograms' field of view 
[S]. Flaring regions are statistically more flux- massive than non- flaring ones. 
As a result, the unsigned flux is considered a "standard" flare forecasting cri- 
terion that, however, is not without limitations (e.g.. iLeka and Barnesl l2003p . 
especially when it comes to flux-massive, but quiescent, active regions. For a 
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Table 1. Values and uncertainties (in parentheses) of the unsigned magnetie flux ^tot, thel 
fractal dimension Do, the inertial-range scaling exponent C(3), and the power-spectrum index 
a for the three nearly simultaneous and coaligned magnetograms of NOAA AR 10930 acquired 
on 11 December 2006 (Figure [3ll. | 



Ref. 


Observation 


Pixel size 


'I'tot 








(Fig.) 


Time (UT) 


(arcsec) 


(xlO^^ Mx) 


Do 


C(3) 


OL 


|3ji 


13:10 - 16:05 


0.158 


2.85 


1.54 (0.04) 


1.35 (0.03) 


3.00 (0.02) 




13:52:01 


0.605 


2.12 


1.43 (0.02) 


1.67 (0.04) 


3.32 (0.04) 




14:27:01 


1.980 


3.60 


1.41 (0.03) 


1.49 (0.03) 


2.24 (0.05) 



parameter to be characterized as having a significant predictive capability, it 
has to perform better than the unsigned flux. 

3. Dependence of Scale-Free and Multi-Scale Parameters on the 
Spatial Resolution 

Prior to analyzing our extensive MDI dataset we perform a sensitivity test of 
the parameters introduced in Section [2.2l to the spatial resolution of the studied 
magnetogram. Ideally, a parameter reliable enough to distinguish flaring from 
non-flaring active regions should be fairly insensitive to varying spatial reso- 
lution. If this were not the case, then one should at least be able to model 
the parameter's variations with changing resolution. If these conditions are not 
fulfilled, likely threshold values of the parameter that one may use to iden- 
tify flaring regions before they flare are resolution-dependent and. as a result, 
instrument-dependent. 

We use the three magnetograms of NOAA AR 10930 shown in Figure |3] to 
test the parameter values on data with varying spatial resolution. The results are 
provided in Table [T] for the parameters $tot [Equation ([5|)], Dq [Equation ([T])], 
Q{q = 3) [Equation ([3])], and a [Equation (|4])]. Our findings can be summarized 
as follows: 

i) Despite coalignment and near simultaneity, the unsigned magnetic flux $tot 
shows distinct differences for the three different magnetograms: the MDI 
low-resolution magnetogram (Figure |3];) shows w 25% larger flux than the 
SOT/SP magnetogram (Figure |3^), while the MDI high- resolution magne- 
togram (Figure [3}d) shows « 28% less flux than the SOT/SP magnetogram. 
Some additional testing, not shown here, has been performed in an attempt 
to explain these large differences. In particular, attempting to degrade the 
highest-resolution magnetograms (SOT/SP and high- resolution MDI) to sim- 
ulate situations of a larger pixel size (2 arcsec and more) we find a very weak 
decreasing trend for the unsigned flux that reaches up to sa 3.5% between 
the highest- (original) and the lowest- (degraded) resolution magnetogram. 
Likely, therefore, the « 30% flux difference is not due to the different spa- 
tial resolution. Similar results arc obtained when the spatial resolution of 
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the SOT/SP and high-resolution MDI data is decreased by resampling. We 
cannot be certain about the source(s) of the discrepancy but, given that we 
have used SP and MDI data of similar calibration levels, one might attribute 
this discrepancy to differences in the observations, processing, and calibration 
between the SP and the two MDI magnetographs. Discussing those differences 
exceeds the scope of this work. 

ii) The values of the scale- free fractal dimension Dq are fairly consistent, despite 
the widely different spatial resolution and the different instruments. One no- 
tices a significant (i.e., beyond error bars) decreasing tendency for decreasing 
resolution but the overall decrease for the « 13-fold difference in resolution 
between the SOT/SP and the low-resolution MDI data is only w 8.5%. 
This is probably because the fractal dimension qualitatively highlights the 
morphological complexity of the studied self-similar structure that is being 
reflected adequately on seeing- free (SOT/SP and MDI) magnetograms largely 
regardless of spatial resolution and magnetic flux content. 

in) The values of the multiscale inertial-range scaling exponent C(3) and the 
power-spectrum scaling index a appear strongly dependent on the spatial 
resolution and/or other instrumental characteristics. The variation of both 
indices is not even monotonic, with their peak values corresponding to the 
intermediate case of the high-resolution MDI magnetogram. One might spec- 
ulate that this happens because different spatial resolution changes quanti- 
tatively, but perhaps not qualitatively, the multiscale character of the data. 
In this sense, despite different parameter values, all three magnetograms of 
NOAA AR 10930 show significant intermittency (C(3) > 1, with lack of 
intermittency reflected on (^(3) = 1) and non-Kraichnan/non-Kolmogorov 
turbulence (a » 3/2 and 5/3, respectively). 

The susceptibility of the multiscale parameters C(3) and a, but not of the 
scale- free parameter Dq, to the spatial resolution implies caution when utilizing 
multiscale parameters to distinguish flaring from non-flaring active regions. At 
the very least, quantitative results in this case should not be generalized to 
different data sets. We therefore stress that the results described in the next 
Sections for the multiscale parameters correspond exclusively to the MDI full- 
disk spatial resolution of « 2 arcsec. 

4. Comparison of Parameters for Flaring and Non-Flaring Active 
Regions 

To determine whether any of the fractal or multifractal parameters Dq, C(3), and 
a, including $tot as a reference, can distinguish flaring from non-flaring active 
regions, we perform two tests: a more stringent one, that compares the preflare 
(96 minutes in advance, at most, per the cadence of the full-disk MDI magne- 
tograms) values of the parameters for flaring active regions to the peak values of 
the parameters for non-flaring regions, and a more liberal one, that compares the 
peak values of the parameters for both flaring and non-flaring regions. Finding 
a distinguishing pattern in the flrst test would mean that the studied parameter 
may have a short-term predictive capability. If the flrst test fails but the second 
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Figure 5. Comparison between tlie preflare values of scale-free and multiscale parameters of 
flaring active regions and the respective pealc values of non-flaring regions. Shown are (a) the 
unsigned magnetic flux ^tot, (b) the fractal dimension Dq, (c) the turbulent power-spectrum 
index a, and (d) the change in the inertial-range scaling exponent C(3). 



test gives some distinguishing patterns, the studied parameter provides clues 
about the expected flare productivity of an active region, without necessarily 
implying when flares will occur. This claim is found in multiple instances in 
the literature {e.g., McAteer, Gallagher, and Ireland, 2005" for Dn : I Abramenkol 



[MJ51 and IHewett et al.[ 120081 for a: IConlon et aU ,2008. for other multifractal 
parameters) . 

4.1. Preflare vs. Peak Non- Flaring Parameter Values 

Figure [5] depicts the preflare values of $tot (Figure [S^), Dq (Figure [S}d), and a 
(Figure [5]:). It also provides the change between the preflare and the postflare 
values of ^(3) (Figure[5ji) for flaring regions, compared with the peak ^(3)-change 
for non-flaring ones. Flaring regions are divided into "M-flaring" (regions that 
have given at least one M-class, but not a X-class, flare; blue squares) and "X- 
flaring" (regions that have given at least one X-class flare; red squares). Since 
peaks of the studied parameter timeseries are taken, each active-region timeseries 
has been inspected separately to remove spurious effects in the parameters' 
evolution. The most prominent source of these effects is contamination due to 
flare emission, in the case of large white-light flares, with temporary instrumental 
or data problems playing a secondary role. In case spurious effects are detected, 
the affected parameter value is replaced by the value interpolated for a given 
time. 

Inspecting Figure [SJ wc first notice that only the ^tot-values of the fiaring 
regions (Figure [5^) show some tendency to occupy the upper <I>tot-rangc in the 



SDLA_ms.tex; 4 January 2011; 1:54; p. 14 



Are Flaring Active Regions More Complex than Others? 



15 



plot. The preflare values of the fractal dimension Dq (Figure [S})) appear to 
be more or less uniformly distributed between « 1.2 and w 1.8, with a mean 
value of w 1.41 and a standard deviation of « 0.08. For non-flaring regions, the 
peak fractal dimension has a mean of w 1.44, with a standard deviation « 0.14. 
Qualitatively, we reproduce the result of|McAteer, Gallagher, and Ireland] (|2005p 



who found Dq > 1.2 for flaring regions, but we also find that Da > 1.2 for all 
active regions in our sample. 

A similar behavior is seen when the scaling index a of the turbulent power 
spectrum is examined (Figure [S};) . Indeed, the prefiare a- values show a mean 
w 2.4 and a standard deviation ?» 0.32, while the peak a-values for the non- 
flaring regions have a mean « 2.60 with a standard deviation w 0.42. No active 
region shows a a- value smaller than the Kolmogorov index of 5/3 and very few 
regions, among them three M-flaring ones, show a < 2. We conclude that a strong 
departure from a Kolmogorov turbulent spectrum is not a characteristic of some 
(flaring) active regions but one of most active regions. Given the a-dependence 
on spatial resolution, however (Section|3]), we cannot be certain about the "true" 
a- value. 

For the multifractal inertial-range scaling exponent ^(3) we have compared 
the change in values between the preflare and the postflare phase with the peak 
change in values for non-flaring active regions (Figure [5ji). This was chosen 



because in Georgoulis (|2005p we inferred a significant, permanent decrease in 
C(3) from preflare to postflare in a single flaring region. Unfortunately, this 
feature does not survive here, where more comprehensive statistics are involved: 
for the 111 major flares included in our sample (24 X-class, 87 M-class) the 
host active regions show a prefiare / postflare increase in C(3) in 59 cases, 
with a decrease between the preflare and postflare C(3)-values inferred in the 
remaining 52 cases. Clearly, the peak ({3) changes do not correspond to fiaring 
active regions with the exception of the two M-flaring regions NCAA AR 9087 
and 10596, that show two of the sharpest C(3)-decreases in the postflare phase, 
with amplitudes 0.66 and 0.31, respectively. If these two regions were studied 
in isolation and one ignored the C-dependence on the spatial resolution, then 
one might have concluded that large flares indeed relate to sharp decreases in 
C(3), meaning a decrease of the degree of the photospheric intermittency in the 
active regions after the flare. Given our spatial-resolution test and the large 
active-region sample, however, such a conclusion is unjustified. 

The apparent failure of the first, stringent test comparing the preflare values 
of the studied parameters for flaring regions with the respective peak values for 
non-flaring ones indicates that none of the scale-free and multiscale parameters 
shows any notable short-term flare prediction capability or, at least, any better 
predictive capability than the conventional unsigned magnetic flux. 

4.2. Peak Flaring vs. Peak Non-Flaring Parameter Values 

We now perform the second test by comparing the peak values of the studied 
parameters for both fiaring and for non-flaring regions. The results of this test 
are shown in Figure El 

Figure [5^. depicts the results for the unsigned magnetic flux $tot • It shows 
more clearly than Figure [S^ that the peak $tot-values tend to occupy the upper 
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Figure 6. Same as Figure [5] but showing the peak, rather than the preflarc, values of the 
examined parameters for the flaring active regions. 



part of the $tot-range. Somewhat more tell-tale is the difference between Figures 
^jp and [GJd in terms of the preflare and the peak fractal dimension, respectively: 
the peak Do-values for flaring regions in Figure also tend to occupy the upper 
part of the Do-range, meaning that flaring regions statistically tend to have 
higher fractal dimension, as |McAteer, Gallagher, and Ireland] (|2005p previously 
reported. 

Regarding the indices a of the turbulent power spectrum (Figure [nt), there is 
no distinct difference with Figure [SJ;. The peak a- values for flaring regions also 
tend to occupy the higher a-range but with larger dispersion than the fractal 
dimension Do. This being said, the active regions with the highest peak a- values 
in our sample happen to be non-flaring ones. 

When the peak (^(3) change for both flaring and non-flaring regions is con- 
sidered (Figure [6ji) we find that flaring regions show a similar pattern with 
non-flaring ones. Only for the X-flaring regions there seems to be a weak sta- 
tistical preference for stronger (^(3)-decreases as compared to increases, with a 
ratio 11:6. The probability that this preference is by chance, however, is rather 
high, of the order 0.09. This probability was calculated by means of a binomial 
probability function, assigning a 0.5 probability that the peak C(3)-change will 
be positive and assuming 11 positive chance "hits" out of 17 independent trials 
(since the ({3) > probability is 0.5 the binomial distribution becomes sym- 
metric with respect to chance hits, so identical results would be reached in case 
we assumed 6 negative chance "hits" out of 17 trials). For M-flaring regions, 
stronger C(3)-decreases barely dominate, with a ratio 32:28. This time, however, 
the binomial probability that this result is random is « 3 x 10"^"*, so we can 
safely conclude that there is no clear preference of C(3)-decreascs over increases 
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Table 2. Summary of moans and standard deviations (in parentheses) for the prefiarc and peakl 
values of ^tot, the scale-free Do, and multiscalc parameters a and C(3), calculated in our sample 
of 17 X-flaring, 60 M-flaring, and 293 non-flaring active regions. The preflare active-region values 
correspond to 24 X-class and 87 M-class flares. The discrepancy between numbers of flaring regions 
and flares is because some flaring regions flare repeatedly over the observing interval. | 





Preflare values 




Peak values 




Active Regions 


X-flaring 


M-flaring 


X-flaring 


M-flaring 


Non-flaring 


Data points 


24 


87 


17 


60 


293 


*tot (xl022 Mx) 


4.57 (1.63) 


3.55 (1.72) 


4.74 (1.68) 


3.84 (1.56) 


2.58 (1.21) 


Do 


1.44 (0.07) 


1.40 (0.08) 


1.51 (0.05) 


1.50 (0.05) 


1.44 (0.13) 


a 


2.38 (0.29) 


2.40 (0.33) 


2.77 (0.36) 


2.73 (0.36) 


2.60 (0.42) 


f (3)-change 


0.012 (0.06) 


-0.001 (0.11) 


-0.03 (0.20) 


-0.04 (0.26) 


-0.004 (0.29) 



in case of M-flaring regions. For non-flaring regions, the peak C(3)-changes are 
almost evenly divided between decreases and increases, with a ratio 146:147, 
that is again not random (the binomial probability that this result is random, 
given the large sample size, is practically zero). Overall, it becomes clear that 
one cannot use the preflare (Figure[S}D) or the peak (FigurejSla) change in C(3) to 
assess the flaring productivity of an active region, meaning that flaring regions 
do not show distinguishable, sharp changes in their degree of intermittency. 

Table [D quantifies Figures [5] and [SI providing means and standard deviations 
for each depicted distribution. It shows that only the peak values of the unsigned 
flux ^tot and the fractal dimension Dq are somewhat different between flaring 
and non-flaring regions but the respective dispersions are such that there is 
considerable mixing of values between the three active-region populations. For 
the preflare values of ^tot and Dq there is even more mixing and, in terms of 
Dq, X- and M-flaring regions are practically indistinguishable. When it comes to 
a- values and changes in C(3), Table [2]- along with Figures [5] and H]- shows that 
these multiscale parameters simply cannot be used to distinguish flaring from 
non-flaring regions, let alone predict large flares within a given time span. 

4.3. Flare Forecasting Probabilities 

In Sections 14. 1114.21 we found that neither the preflare nor the peak values of our 
scale-free and multiscale parameters seem capable of distinguishing flaring from 
non-flaring regions. The unsigned flux $tot tends to score better than the scale- 
invariant Do, the multiscale a, and the change in C(3). To further quantify these 
results we calculate here the conditional probability of an active region being 
a flaring one if a given parameter inferred from any one of its magnetograms 
exceeds a preset threshold. A way to do this is by using Laplace's rule of suc- 
cession and proceeding to a Bayesian inference of the predictive probability as 
follows: assume that F magnetograms of flaring regions and A'^ magnetograms 
of non-flaring ones exhibit a value i? of a parameter that exceeds a threshold 
^thrcs- Then, the conditional probability p that an active-region with R > i?thros 
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Figure 7. Conditional probabilities of active regions to liost a major flare, cither X-class (left 
column) or M-class (right column), with respect to the normalized (against the maximum) 
threshold of a given parameter, with ^tot (blue squares), Dq (red squares), and a (green 
squares) examined. The upper (a, d) and middle (b, e) rows provide the 12- and 24-hour 
conditional probabilities, respectively. The lower row (c, f) provides the probability without 
a time limit, that is, the probability of major flaring at a future time when the active region 
is still visible in the disk. The maximum values against which the thresholds were normalized 
are 8.8 X 10^^ Mx, 1.65, and 4.17 for $tot, A), and a, respectively. 



will be a flaring one is given by ( |Jaynes[ 120031 pp. 155-156) 



F+l 
N + 2 



with uncertainty 5p 



N + 3 



(6) 



This probability rule was also used by IWheatlandl (I2005P to test another solar- 
flare prediction method. To compare directly between the conditional probabil- 
ities of the various parameters, we normalize the thresholds with respect to the 
maximum value of each parameter appearing in Figures [5] and |6l We examine 
only three of the four parameters - <I>tot, Dq, and a - because the changes of C(3) 
obviously show very similar patterns for flaring and non-flaring regions. Using 
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the NOAA/GOES and Yohkoh/RXT flare catalogs (Section 2.1), we consider a 
given active-region magnetogram as a preflare one if a major flare happened in 
the region within a given, preset timeframe from the time that the magnetogram 
was recorded. Figm^e [7] shows conditional probabilities for two timeframes: 12 
hours (Figure [7^, d) and 24 hours (Figure [7}d, e), while Figures [71; and f show 
conditional probabilities without any timeframe assigned: if an active region 
flared at any time after the magnetogram was taken, then this magnetogram 
is considered a preflare one. Notice that the flares under examination may have 
occurred when a given region has moved beyond the 30° E-W meridional zone of 
analysis. This does not affect the analysis, however, as the preflare magnetograms 
were acquired when the region was still within the analysis zone. 

All plots of Figure [7] illustrate that the unsigned flux ^tot, a conventional 
activity predictor, is generally more effective in predicting major flares than both 
the scale-free fractal dimension Dq and the multiscale turbulent power-spectrum 
index a. Differences between $tot and (Dq, a) are smaller (but also reflect 
small flare probabilities, of limited practical use) in case of the most demanding 
prediction, the one with a 12-hour timeframe for flares >X1.0 (Figure[7^). In all 
other cases $tot gives much higher (well beyond error bars) probabilities than Dq 
and a. The predictive ability of Dq appears comparable with, or slightly higher 
than, that of $tot only for normalized thresholds -Rthros > 0.9 for the 12- and 
24-hour prediction timeframes. For the same timeframes, the predictive ability 
of $fot drops for -Rthros 0.8 (Figures [7^, b) and i?thrcs ^ 0.6 (Figures [TJi, e). 
This is because the upper $tot-ranges in these cases are occupied by non-flaring 
(within the preset timeframes) active regions - this is one of the limitations 
for using the unsigned flux as a flare predictor. The power-spectrum index a 
exhibits similar behavior, but for higher i?tiircs ^ 0.8, in all plots of Figure [71 
Comparing the multiscale a with the scale- free Dq, we note that a works some- 
what better, especially for larger timeframes. This is in line with lAbramenkol s 
(|2005p suggestion that a better reflects future flare productivity. However, recall 
that a depends sensitively on the spatial resolution of the observing instrument, 
contrary to Dq (Table [1]) . Hence the results of Figure [71 concerning a should be 
viewed as holding exclusively for MDI low-resolution magnetograms. 

5. Summary and Conclusion 

This study investigates previous claims on the efficiency of fractal and multi- 
fractal techniques as reliable predictors of major solar flares and/or parameters 
reflecting the overall flare productivity of solar active regions before they actually 
flare. From the array of parameters implemented in the literature, we select three 
of the reported most promising ones: the fractal dimension, the multifractal 
intermittency index, and the scaling index of the turbulent power spectrum. 
Our objective is not to judge the methods per se but, rather, to test the notion 
of utilizing fractality and multifractality to gain predictive insight into major 
solar flares. 

Statistical analyses such as this one must guarantee that the assembled active- 
region sample is representative: the sample must contain numerous flaring and 
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non-flaring regions. Comprehensive statistics often help avoid the interpretation 
of incidental signals as statistically significant behavior. Section HTTI (Figure ISH) 
includes examples of results that might have been interpreted in a misleading 
way had the statistics of our active-region sample been insufficient. 

We study 370 SOHO/MDI low-resolution (1.98" per pixel) timeserics of active-| 
region magnetograms, 293 of which correspond to active regions without major 
flares and 77 correspond to M- and X-class flaring regions. MDI line-of-sight fields 
are used for regions within 30° of the central meridian in order to approximate 
the longitudinal-field component with the normal-field component and avoid any 
corrections or otherwise modifications of the original MDI data. We find that 
neither scale-free (fractal) nor multiscale (multifractal) techniques can be used 
to predict major fiares, or for the a priori assessment of the fiaring productivity 
of active regions. In particular, wc find that their diagnostic capability is not 
better than that of the unsigned magnetic flux of active regions, a traditional, 
but unreliable, activity predictor. Since the fractal and multifractal measures 
tested here are less effective than the unsigned flux (Figure [T]), they should not 
be used for flare prediction or for flaring productivity assessment. 

On the fundamental question of whether flaring active regions are more frac- 
tal, multifractal, or turbulent than other, non-flaring ones, the answer per our 
results has to be negative: flaring regions tend to exhibit relatively large peak 
values of scale-free and multiscale parameters but these values, or even higher 
ones sometimes, arc also exhibited by non-flaring regions. For all statistical 
distributions, the means and standard deviations are such that the different 
populations of flaring and non-flaring regions overlap considerably (Tabic [2]). 

At this point we emphasize our willingness to follow the guidelines of multiple 
previous studies in the inference of the above fractal and multifractal parameters. 
In particular, we followed'McAtcer, Gallag her, and Ireland] (j2005p when inferring 
the fractal dimension Do , Abramenko (2005]) when inferring the turbulent scaling 
index a (despite the fact that Abramenko worked exclusively on high-resolution 
MDI magnetograms), and a previous work of this author ( [Georgoulis, 2005] ), 
together with lAbramenko et aTl (j2003p , when inferring the intermittency index 
C(g). As a result, the findings of both [McAteer, Gallagher, and Ireland (I2005P 
and lAbramenko] (j2005p were qualitatively reproduced in this analysis, while we 
showed that the distinct C(3)-behavior reported by |GeorgouIis] (|2005p was just 
one incidental case and not part of a systematic trend. 

In addition, this work (Section [3]) exposes a dependence of multiscale pa- 
rameters Cisi) and a on the spatial resolution of the studied magnetograms. In 
contrast, the scale-free Dq appears fairly insensitive to varying spatial resolution. 
Therefore, results and comparisons for ^(3) and a in Section 2] are valid only 
for MDI low-resolution data and should not be generalized to data sets of other 
instruments. Possible susceptibility of the Dq- value should also be studied with 
respect to the threshold it requires, unlike ^(3) and a. This investigation has not 
been carried out here. In previous works, however. IMeunierl(|1999p reported a de- 
creasing trend of Dq with increasing threshold, while [Janfien, Vogler, and Kneer| 
(|2003p reported a slighter decrease, or a near insensitivity, of for increasing 
thresholds, in case these thresholds arc sufficiently above noise levels or the 
magnetic field data have been treated for noise, respectively. 
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It is useful to mention here a very recent result by |Abramenko and Yurchyshyn[ 
(|2010p that the turbulent power-spectrum index a, either alone or coupled with 
the integral of the power-spectrum for all wavenumbers, correlates better than 
*I*tot with the flaring index in a large sample of 217 active regions recorded 
in high-resolution MDI magnetograms. While correlating some parameter with 
the flaring index is not identical to inferring the predictive capability of this 
parameter, these results appear in likely contrast with the results presented here. 
Further investigation is clearly needed, therefore. Nonetheless, some convergence 
of views appears in that multiscale parameters may not be ideal tools for solar 
flare prediction (Abramenko, 2010, private communication). 

Perhaps more instructive than pointing out the inability of scale-free and 
multiscale techniques to assess a priori the flaring record of active regions is to 
explain why this is the case. In this author's view, there are at least two distinct 
reasons that justify our findings: 

First, fractality and multifractality are extremely widespread in the solar 
atmosphere, eruptive and quiescent alike. This may well be due to the turbulence 
dominating the magnetic-flux generation and emergence process (see Introduc- 
tion). For example, recall the fractality of white-light granules ( [Roudier and Muller, 1987H 
|Hirzberger et ai, 1997[ ), the fractality and multifractality of active regions and 
the quiet- Sun magnetic field (Schrijver et ai, 1992||Cadavid et al., 1994 [Meunier, 1999H 



|Janfien , Vogler, and Kneer, 2003), the fractality of flares and sub-flares in the 
EUV (jAschwAndon Parncll, 2002; Aschwanden and Aschwandeiill2008all2008bp .| 



the fractality of the quiet network in the EUV (Gallagher et ai, 1998[), that 



of Ellerman bombs in off-band Ha (Georgoulis et ai, 2002 1, and others. The 



fractal dimension in most, if not all, of these works varies between 1.4 and 
1.8, practically indistinguishable from the fractal dimension of active regions 
found here. As a result, it appears unlikely that these same methods may reflect 
particular characteristics of active regions, let alone flare productivity. 

Second, there is a lack of correlations between the fractal dimension in the 
photosphere and that of the overlying chromosphere and corona, where major 
flares occur. [Dimitropoulou et al.\ (j2009p assumed nonlinear force-free magnetic 
fields extending above the photosphere and calculated volumes of enhanced 
electric currents and steep magnetic gradients from these extrapolated fields. 
They found no correlation between the three-dimensional fractal dimension of 
these volumes and that of the two-dimensional photospheric boundary. In other 
words, all photospheric "memory", in terms of fractality and multifractality, 
is erased above the photosphere due to the fact that these unstable volumes 
become nearly space-filling slightly above this boundary. Attempting to assess 
the fractality of layers higher than the photosphere - where fiares occur - by 
using the photospheric fractality as a proxy will not yield meaningful results, 
similarly to the lack of correlation between photospheric electric currents and 
the coronal X-ray brightness ( [Metcalf et ai, 1994[ ). 

In addition, it is possible that both flaring and non-fiaring regions share 
a similar degree of self-organization in the distribution of their magnetic free 
energy, as reported by [Vlahos and Georgoulis (|2004p . Flaring regions have an 



"opportunity" to show their self-organization via fiaring, with fiares inheriting 
the statistics of their host active regions, while non-flaring regions retain this 
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property without demonstrating it. In this sense i) fractaUty alone cannot be 
responsible for flaring, and ii) fractality, as a global characteristic of the active- 
region atmosphere, cannot be used to determine a priori which active regions 
will flare. 

There are, of course, sophisticated multiscale techniques not treated in this 
work, such as wavelet methods used to extract the magnetic-energy spectrum in 



active regions (Hewett et al., 2008) or to distinguish active regions from quiet 



Sun for further treatment (Conlon et al., 2010), or the flatness function and 



its intcrmittency index (Abramenko, Yurchyshyn, and Wang, 2008). While we 



cannot comment on methods that we have not tested, per our conclusions it 
would seem rather surprising if a scale-free or multiscale technique delivered a 
notable improvement in our forecasting ability, as this would apparently contra- 
dict what scale-free and multiscale behavior caused by self-organization is meant 
to imply: spontaneity in the system's dynamical response to external forcing, 
both in timing and in amplitude, and hence a lack of certainty in predicting this 
response. 

Let us finally mention that alternative flare prediction approaches have been 
developed in recent years. Rather than fractality, multifractality, or intcrmit- 
tency and turbulence, these methods rely on parameters stemming from mor- 
phological and topological characteristics of active regions, such as those of 



general (Falconer, Moore, and Gary| 120061 |Schrijver 120071 


Georgoulis and Rustl 


120071 ILeka and iiarnesl 120071 IJVlason and Hoeksemal 12010 


), or those of the sub- 


surface kinetic helicity prior to active-region emergence ( 


Reinard et al, 2010), 



advances in the forecasting of major solar eruptions or whether forecasting will 
remain inherently probabilistic which, per our results, seems entirely possible. 
In any case, fractal and multifractal methods - perhaps not extremely use- 
ful as eruption predictors - will always be excellent tools for a fundamental 
understanding of the origins and nature of solar magnetism. 
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